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Abstract 

We propose new methodologies for stabilizing all-metal antiaromatic clus- 
ters like: Al 4 Li 4 . We demonstrate that these all-metal species can be stabi- 
lized by complexation with 3d-transition metals very similar to its organic 
counterpart, C4H4. Complexation to transition metal ions reduce the fron- 
tier orbital energies and introduces aromaticity. We consider a series of such 
complexes [?7 4 (Al4Li4)-Fe(CO)3, r/V^ALiL^-Ni and (AULi^Ni] and make a 
comparison between the all-metal species and the organometallic compounds 
to prove conclusively our theory. Fragmentation energy analysis as well as 
NICS support similar mechanism of complexation induced stability in these 
all-metal molecules. 
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The concept of aromaticity and antiaromaticity is of fundamental im- 
portance in chemistry. From the simple Hiickel theory to the more refined 
concepts like ring-currents and critical point charge- densities, the field has 
evolved over a period of half-a-centuryPQ. The idea has been extended from 
initially a small class of organic 7r-conjugated systems to now many inorganic 
molecules and molecular clusters[2], with rapid experimental verifications 
through actual synthesis and characterizations [3J. 

The recent report of the first all-metal antiaromatic complex, Al 4 Li 4 
shows the generalizations and usefulness of this idea[3]. Also, in the last few 
years there have been reports of aromaticity in all-metal clusters^- However, 
unlike their organic counterpart, C4H4, where the energy separation between 
the a and n orbitals is substantial, these all-metal molecules have closely 
placed orbitals and thus have poor a-n separation. Mostly, the aromatic 
characteristic is associated with delocalized electrons (n electrons). However 
recently it has been reported that these clusters are more a-aromatic than 
7r-antiaromaticjn]. As a result, there has been a confusion whether to call 
these complexes aromatic or antiaromatic. The controversy can only be set- 
tled after the successful synthesis followed by unambiguous crystal structure 
determinations. For a stable molecular crystal, measurement of bond lengths 
alternations as well as the charge densities (ring critical points) are well es- 
tablished parameters for characterizing aromaticity /antiaromaticity [7j. 

These clusters however have till now been synthesized only in the gas- 
phase by laser vaporization technique and is thus insufficient in providing 
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structural details. Recently, we have shown that these materials, A1 4 M 4 
(M=Li, Na and K) are also very good candidates for higher-order nonlin- 
ear optical (NLO) applications due to charge transfer from the highly elec- 
tropositive ion (Li) to the Al 4 -ringsjS]. Unfortunately, the stability is one 
issue which hinders any applications let alone firm establishment of the basic 
understanding. 

The synthesis of antiaromatic molecules is difficult because of their insta- 
bilities. Cyclobutadiene, (C4H4), a An electron system remained non-isolated 
for a longtime before Longuet-Higgins and Orgel proposed in a landmark 
paper, the concept of stabilization through complexation with a transition 
metal to form an organometallic compoundjO]. The compound was synthe- 
sized soon after [TO]. In the following, we justify this simplistic model for the 
small Al 4 -clusters and propose a few very stable complexes of these all-metal 
species. Parallelly, we also compare and contrast the energetics with their 
organic analogues (C 4 H 4 complexes). 

We have performed a closed shell calculation for singlet and an open shell 
calculation for the triplet state at the 6-311G(d,p) basis set level. Elec- 
tron correlation has been included according to the DFT method using 
Becke's three parameter hybrid formalism and the Lee- Yang-Parr function- 
al (B3LYP) available in the Gaussian electronic structure set of codes[TT|. 
The geometries obtained from the B3LYP method have been shown to be in 
very good agreement with the measured photoelectron spectra in such small 
clusters |T2]. 
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In Fig. 1, the minimum energy ground state structures are shown for 
C4H4 and Al4Li4. Simple Hiickel 7r-electron theory predicts a triplet square 
geometry for C4H4 with equal C-C bond lengths However inclusion of 
interaction with the underlying a backbone stabilizes the C4H4 molecule 
in a singlet state with rectangular geometry. This is a good example of 
Jahn- Teller distortion or Pierls instability in low-dimensional system which 
allows stabilizations through bond length alternation. In fact, in this picture, 
the square geometry actually corresponds to a transition state between two 
degenerate rectangular ground state structures. In Table 1, the total energies, 
bond length alternation (Ar, defined as the average difference between the 
bond lengths of two consecutive bonds in the 4-membered ring) for both the 
states are tabulated. The rectangular C4H4 [A(i)] is more stable than the 
square geometry [A(ii)] by 6.2 Kcal/mol. Thus, a triplet square geometry is 
expected to be the transition state for processes such as ring whizzing, where 
one rectangular form is converted into the other (an in-plane rotation of 90°) , 
in harmony with time-resolved transition state studies for the tub-inversion 
in l,3,5,7-cyclooctatetraenejH|. 

Since in Al 4 Li 4 the a-n separation is poor, the Hiickel 7r-electron picture is 
completely invalid. In fact, the n electrons in this case interacts more strongly 
with the a backbone and we expect a distorted structure as the ground state. 
The ground state structure for the singlet state is found to be a rectangular 
Al4-geometry with surrounding Li atoms, forming a C2/1 symmetry group 
[B(i)]. The same structure has been found in previous calculations as welljH]. 
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Another low-energy structure for Al 4 Li 4 found by optimizing the geometry 
for the singlet square geometry is a diamond shaped structure. It has a 
D 2 h symmetry [B(ii)] and is 15 Kcal/mol higher in energy than the stable 
C 2 h geometry [B(i)]. Thus, the rectangular Al 4 -ring in [B(i)] corresponds 
to the ground-state geometry for Al 4 Li 4 . The structural distortion leading 
to a magnetic triplet state with D ih symmetry [B(iii)] is found to lie 55 
Kcal/mole above the ground state singlet [B(i)]. There also exist a low- 
energy structure with the same geometry as the ground state (C 2 /i symmetry) 
with purely spin excitation at an energy only 5 Kcal/mole above the ground 
state. This triplet geometry [B(iv)] for Al 4 Li 4 does not have a counterpart 
in C 4 H 4 , clearly explaining the low a-n gap and the existance of continuum 
like metallic states in Al 4 Li 4 . 

Existence of a very stable rectangular ground state structure together 
with a square geometry as the transition state for the Al 4 -ring similar to those 
for C 4 H 4 suggests that these inorganic clusters are anti-aromatic in nature. 
This is expected since, the highly electropositive Li atoms donate electrons to 
the Al-atoms, thereby creating a species of the type Al 4 4_ , isoelectronic with 
C 4 H 4 . Thus we can safely consider Al 4 Li 4 as a 4n electron system with the 
7T-HOMO (highest occupied molecular orbital) being a non-bonding molec- 
ular orbital as like in C 4 H 4 . In the context of C 4 H 4 , Longuet-Higgins sug- 
gested that such a system can be stabilized if the non-bonding electrons form 
bonding molecular orbitals with suitable low energy d-orbitals of a transition 
metal. For this to happen however, the energies of the d-orbitals should lie 
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close to the low-energy levels of the molecule alone. In the following, we pro- 
pose a few stable complexes of Ai4Li4 and compare their formation energies 
in comparison with their organic analogues. 

1 Fe(CO)3 complex: 

A molecular complex, 77 4 (C4H 4 )-Fe(CO)3, has been recognized through the 
formation of such bonding molecular orbitals and this complex has been 
reported to be quite stable [To]. In fact, oxidation of this complex releases the 
C4H4 ligand which is a stable source for the highly reactive cyclobutadiene 
in organic synthesis [IS]. For the Ai4Li4, we perform ground state energy 
analysis on the similar system, 7] 4 (Al 4 Li 4 )- Fe(CO)3, using the same level of 
theory mentioned above. Both 77 4 (Al 4 Li 4 )-Fe(CO)3 and its organic analogue 
have substantial stability (see Fig. 2 for structures). Al 4 Li 4 indeed forms a 
stable jf complex with Fe(CO)3 [2(h)]. The stability of the complexes are 
investigated using the following fragmentation scheme: 

7? 4 (C 4 H 4 )-Fe(CO)3 = C4H4 + Fe(CO) 3 

r/ 4 (Al 4 Li 4 )-Fe(CO)3 = Al 4 Li 4 + Fe(CO) 3 

The binding energy for 77 4 (Al 4 Li4)-Fe(CO)3 is found to be 106.04 Kcal/mol 
while that for ^ 4 (C4H 4 )- Fe(CO)3 it is 78.44 Kcal/mol. The comparable bind- 
ing energies for the two compounds suggest that ALil^ is very well stabilized 
in the complex, in-fact even more stabilized than C4H4. 
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The HOMO of the molecule (a non-bonding MO for C4H4), interacts with 
the low-energy d-orbital of the ligand to form a bonding combination in the 
complex. Thus, C4H4 that initially possessed 47r electrons now has two more 
electrons, forming a species of the type C4H 4 2_ , an aromatic molecule. Simi- 
larly, for AI4LI4, the complexation converts it into AUL^ 2- , a well-established 
aromatic complex[S]. The HOMO energies for C4H4 in the free and in the 
coordinated form (derived by performing a single-point energy calculation on 
the C4H4 fragment in the optimized complex) are —0.198 au and —0.157 au 
respectively, while the same for the complex, r/ 4 (C 4 H 4 )-Fe(CO)3, is —0.250 
au. Similarly, for AI4L14, the free and the coordinated form have HOMO ener- 
gies at —0.128 and —0.104 respectively and the complex 77 4 (Ai4Li4)-Fe(CO)3 
has the same at —0.168 au. The stabilization of the frontier orbitals in the 
metal complex in both the systems has its manifestation at the formation of 
the stable structure. The similarity in the difference between the HOMO en- 
ergies of free and coordinated structures and with that of the corresponding 
complex for each molecule suggest that a similar mechanism is operative in 
lowering the energies of the frontier orbitals in stabilizing the complexes. 

The above explanation is based on a simplistic picture of the interaction. 
To verify that indeed such a scheme is valid for a molecule with poor a-n 
separation like Al4Li4, we compare the Ar for both C4H4 and AI4LI4 in the 
free geometry and when they are complexed with the transition metal. For 
C4H4, the Ar is 0.24 A in the free state. In the complex, 7? 4 (C4H4)-Fe(CO)3, 
the Ar for the C4H4 ring is only 0.005 A. Thus, C4H4 when complexed is a 
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square rather than a rectangle and as expected from the 7r-only interaction, 
it behaves as aromatic C4H4 2- . For Al4Li4, Ar=0.13 A in the free state 
while in the complex it is only 0.03 A. This clearly supports that the Al 4 Li 4 
has been converted into Al 4 Li 4 2_ , accounting for its substantial stability due 
to aromaticity. The complexation induced metalloaromaticity in AI4LI4 is 
schematically shown in Fig. 3. While a square (triplet) Al 4 Li 4 is much 
higher in energy than the rectangular Al 4 Li 4 , this square form is stabilized 
on complexation to a transition metal. Same is the case for C4H4, where 
the square form becomes stabilized upon complexation. This is similar to 
the origin of aromaticity in benzene, where the 7r-delocalized structure 
corresponds to a energy minima between two bond-altered Kekule forms with 
D^h symmetry. 

An even more clear picture is derived by performing a calculation for 
the nucleus-independent chemical shift (NICS) [T7j at the GIAO-B3LYP /6- 
311+G(d,p) level. We calculate the NICS at the center of the AI4 ring before 
and after complexation with the Fe(CO)3. For comparison, the same values 
are also calculated for C4H4. In C4H4, NICS values before (C4H4) and after 
complexation (C4H 4 2 ~) are 23.55 ppm and -15.37 ppm respectively. The 
change in sign clearly shows the transition from antiaromatic to aromatic 
nature upon complexation. For AI4LI4, the NICS values change from -11.01 
ppm in the free state to -25.44 ppm on complexation in ^(AULi^)- Fe(CO)3. 
The initial negative magnitude for NICS in free Al4Li4 supports the claim by 
Schleyer et al that ALiLi4 has higher u-aromaticity than 7r-antiaromaticity[6 j . 
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But, the increase in NICS value with same negative sign suggests increased 
aromaticity in these clusters upon complexation, which is expected from the 
7r only picture of the conversion of Al 4 Li 4 to Al 4 Li 4 2 ~. Thus, complexation 
with Fe(CO)3 induces metalloaromaticity in Al 4 Li 4 and thereby stabilizes the 
complex, 77 4 (Al 4 Li 4 )-Fe(CO) 3 . 

2 bis nickel(II) chloride complex: 

Next we consider another very well known example of a stable C 4 H 4 com- 
plex, bis(cyclobutadiene nickel(II) chloride). The tetramethyl derivative for 
the complex crystallizes in a P21/c point group and has a good resolution 
(R=7.0%), CCD reference code, NCBNIB [TS]. We have obtained the struc- 
ture from the database and the methyl groups was substituted by H for 
easy comparison with the Al 4 Li 4 derivative, bis(Al 4 Li 4 nickel (II) chloride). 
Both the structures, bis(cyclobutadiene nickel(II) chloride) and bis(Al 4 Li 4 
nickel (II) chloride) were optimized at the same level of theory as mentioned 
above. The structure for bis(cyclobutadiene nickel(II) chloride) remains sim- 
ilar to that found from the crystal structure. Fig. 4 shows the structures for 
the two complexes. For the organometallic complex 4(i), the bond length al- 
ternation in the C 4 H 4 ring is only 0.05 A. Therefore this bridged chlorine sys- 
tem also shows strong mixing of the d-orbitals from Ni and the non-bonding 
electrons of C 4 H 4 . 

For the all-metal complex 4(h) however, the 7? 4 (Al 4 Li 4 )-Ni binding mode 
is converted into ?7 2 <7 2 (Al 4 Li 4 )-Ni upon optimization. In Fig.5, we show this 
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change in the binding mode. Two of the initial 7r-bonds between Al and 
Ni (Al-Ni distance of the order 2.42 A) are now converted into strong a 
bonds (Al-Ni distance of the order 2.25 A) in the optimized geometry. Such 
a change in bonding pattern from 7r character to a character is quite well 
known in organometallic chemistry [HI]. It is interesting to note that a similar 
phenomenon occurs in the all metal complexes as well. 

A fragmentation analysis on these two molecules is given below. 

bis(cyclobutadiene nickel(II) chloride) = 2 C4H4 + bis(nickel(II) chloride) 

bis(Al 4 Li 4 nickel(II) chloride) = 2 Al 4 Li 4 + bis(nickel(II) chloride) 

The stabilization energies are 68.39 Kcal/mol for 4(i) and 286.77 Kcal/mol 
for 4(ii). Such high stability in the Al 4 Li 4 complex is due to the formation 
of two strong Al-Ni a bonds as discussed above. 

3 Metal sandwich complex: 

Another well known methodology in stabilizing a molecule is to form a sand- 
wich type of geometry where two molecular species can share interaction with 
a transition metal: cyclopentadiene is stabilized in such a geometry resulting 
in the Ferrocene structure [20]. For C 4 H 4 , a simple effective electron number 
(EAN) counting shows that the metal in between the two ligands should have 
10 valence electrons in stabilizing a sandwich of the type:(C 4 H 4 )2M. The sim- 
plest metal with 10 electrons in the valence shell is Nickel(O). Elements in 
the same group like Pd or Pt have a strong spin-orbit coupling and prefer 
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square-planar geometry (16 electron geometry). Thus a coordination number 
of 8 as required in a sandwich complex is not possible with Pd or Pt. After 
performing the geometry optimization at the same level of theory discussed 
above, we find that the structure for (C4H 4 ) 2 Ni is indeed a sandwich geom- 
etry with the two C4H4 rings above and below the Ni atom (see Fig. 6 (i)). 
In this complex, the Ni atom sits symmetrically inside the cavity of the two 
C4H4 rings with a distance of 1.99A from each C4H4 ring. The two C4H4 are 
staggered to each other. 

Similarly, we have been able to stabilize the Al 4 Li 4 cluster by introducing 
it in a sandwich of the type: (AL^Li^Ni. The geometry is shown in Fig. 6 
(ii) (the optimization is performed at the B3LYP/LANL2MB level followed 
by energy calculation at B3LYP/6-311G(d,p) level). The central Ni atom 
sits unsymmetrically in the cavity of the two Al 4 Li 4 rings. A very recent the- 
oretical study on its aromatic analogue, Al 4 2 ~, support our calculations fIT\. 
Interestingly, the Al atoms in the rings bend towards the Ni atom and the 
planarity of the AI4 ring is thereby lost. This is understood from the fact 
that when the 4-7T electrons of each of the two Al 4 Li 4 rings interact with the 
central Ni atom, the requirement of the Al atoms to be in plane with the Li 
atom is no longer important. Instead, the sandwich like structure with 18 
electrons gives an extra stabilization keeping the whole system electrically 
neutral. The stability of these complexes are investigated using the following 
fragmentation scheme: 

(C 4 H 4 ) 2 Ni = 2 C 4 H 4 + Ni(0) 
11 



(Al 4 Li 4 ) 2 Ni = 2 Al 4 Li 4 + Ni(0) 

where Ni(0) is in a 3 F state. Al 4 Li 4 binds strongly to the Ni(0) and has a 
binding energy of 134.285 Kcal/mol. For C 4 H 4 this binding energy is 150.819 
Kcal/mol. 

In conclusion, we have demonstrated for the first time that the all-metal 
species like Al 4 Li 4 can be stabilized by complexation with 3d-transition met- 
als, very similar to its organic counterpart, C 4 H 4 . Although such a complex- 
ation induced metalloaromaticity is a well established concept in the realm 
of organometallic chemistry, we have demonstrated that it is a very general 
and elegant concept which can be used for all metallic molecules with or- 
bitals that are close in energy with the d-orbitals of the transition metal. We 
have also shown that these all-metal complexes have similar binding ener- 
gies and properties like their organometallic counterparts and thus should be 
considered as very good candidates for experimental synthesis. Also, such a 
stabilization will provide a very precise answer to the question of aromatic- 
ity/antiaromaticity and will lead to novel applications of these clusters. We 
believe that our work will motivate synthesis of these molecules as were the 
case in the previous generation for the organometallic complexes. 
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Table 1: The total energies (in au) and the bond-length alternation, Ar (in 
A) for C4H4 and Al 4 Li 4 in different spin states corresponding to different 
lowest energy structures. 



Molecule 


Symmetry 


Spin-state 


Energy 


Ar 


C 4 H 4 (A(i)) 




Singlet 


-154.718 


0.240 


C 4 H 4 (A(ii)) 


D 4ft 


Triplet 


-154.708 


0.000 


Al 4 Li 4 (B(i)) 


C211 


Singlet 


-999.932 


0.130 


Al 4 Li 4 (B(ii)) 


D 2A 


Singlet 


-999.908 


0.120 


Al 4 Li 4 (B(iii)) 


D 4A 


Triplet 


-999.844 


0.000 


Al 4 Li 4 (B(iv)) 


C2h 


Triplet 


-999.926 


0.200 
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Figure 1: Equilibrium minimum energy geometries for C4H4 and A^L^ in 
the singlet and triplet states (See Table. 1). Bond lengths are given in A. 
Ball color: Black=C, White=H, Pink=Li, Light orange=Al. 




Figure 2: Equilibrium minimum energy geometries for (i) ?7 4 (C4H4) -Fe(C0)3 
and (ii) ?] 4 (Al4Li4)-Fe(CO)3. Bond lengths are in A. Ball color: Red=0, 
Violet=Fe. 
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Figure 3: Schematic representation of (i) change in geometry for ring 
whizzing and complexation to transition metal center for A=A1 in Al4Li4 
(A£ 1= 55 Kcal/mol, A£ 2 =100 Kcal/mol); A=C in C 4 H 4 (A£ 1= 6.2 
Kcal/mol, AE 2 —78A Kcal/mol) (ii) Ring whizzing in benzene. 
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Figure 4: Equilibrium minimum energ^ geometries for (i) bis(cyclobutadiene 
nickel(II) chloride)and (ii) bis(Al 4 Li 4 nickel(II) chloride). Distances are in A. 
Ball color: Green=Cl, Blue=Ni. 
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Figure 5: Conversion from rf binding mode for AI4LI4 to rf - cr 2 in the complex 
in Fig. 4(ii). 
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E=- 18 17 878 au 
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Figure 6: Equilibrium minimum energy geometries for (i) (C^^Ni and (ii) 
(Al4Li4)2Ni. Distances are in A. 
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